(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 11.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[      6656,        176]
NotebookOptionsPosition[      6184,        155]
NotebookOutlinePosition[      6538,        171]
CellTagsIndexPosition[      6495,        168]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell["S(k) integral evaluation to calculate theory curves :", "Text",
 CellChangeTimes->{{3.839727351791481*^9, 3.839727364221589*^9}, {
  3.839735477559823*^9, 3.839735523349258*^9}, {3.839735760537407*^9, 
  3.839735767652286*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"fugacity", "[", "tovertf_", "]"}], ":=", 
   RowBox[{"z", "/.", 
    RowBox[{"FindRoot", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"PolyLog", "[", 
        RowBox[{"3", ",", 
         RowBox[{"-", "z"}]}], "]"}], "+", 
       RowBox[{"1", "/", 
        RowBox[{"(", 
         RowBox[{"6", " ", 
          RowBox[{
           RowBox[{"(", "tovertf", ")"}], "^", "3"}]}], ")"}]}]}], ",", 
      RowBox[{"{", 
       RowBox[{"z", ",", "1000"}], "}"}]}], "]"}]}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"kscal", "[", 
    RowBox[{"tovertf_", ",", "efoverer_"}], "]"}], ":=", 
   RowBox[{"2", " ", 
    RowBox[{"Sin", "[", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"Pi", " ", 
        RowBox[{"24", 
         RowBox[{"(*", 
          RowBox[{"observation", " ", "angle"}], "*)"}], "/", "180"}]}], 
       ")"}], "/", "2"}], "]"}], " ", 
    RowBox[{"Sqrt", "[", 
     RowBox[{"1", "/", "efoverer"}], "]"}], " ", 
    RowBox[{"Sqrt", "[", 
     RowBox[{"1", "/", "tovertf"}], "]"}]}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"sls", "[", 
    RowBox[{"tovertf_", ",", "efoverer_"}], "]"}], ":=", 
   RowBox[{"1", "-", 
    RowBox[{"8", " ", 
     RowBox[{
      RowBox[{"tovertf", "^", "3"}], "/", "Pi"}], " ", 
     RowBox[{"NIntegrate", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"a", "^", 
         RowBox[{"(", 
          RowBox[{"3", "/", "2"}], ")"}]}], " ", 
        RowBox[{"1", "/", 
         RowBox[{"(", 
          RowBox[{"1", "+", 
           RowBox[{
            RowBox[{"1", "/", 
             RowBox[{"fugacity", "[", "tovertf", "]"}]}], " ", 
            RowBox[{"Exp", "[", 
             RowBox[{"a", "+", 
              RowBox[{"y", "^", "2"}]}], "]"}]}]}], ")"}]}], "*", 
        RowBox[{"1", "/", 
         RowBox[{"(", 
          RowBox[{"1", "+", 
           RowBox[{
            RowBox[{"1", "/", 
             RowBox[{"fugacity", "[", "tovertf", "]"}]}], " ", 
            RowBox[{"Exp", "[", 
             RowBox[{"a", "+", 
              RowBox[{
               RowBox[{"(", 
                RowBox[{"y", "+", 
                 RowBox[{"kscal", "[", 
                  RowBox[{"tovertf", ",", "efoverer"}], "]"}]}], ")"}], "^", 
               "2"}]}], "]"}]}]}], ")"}]}]}], ",", 
       RowBox[{"{", 
        RowBox[{"y", ",", 
         RowBox[{"-", "Infinity"}], ",", "Infinity"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"a", ",", "0", ",", "Infinity"}], "}"}]}], "]"}]}]}]}], 
  ";"}]}], "Input",
 CellChangeTimes->{{3.804740235757255*^9, 3.804740271980023*^9}, 
   3.816707984185689*^9, {3.817160230777093*^9, 3.817160231467679*^9}, {
   3.817306117280761*^9, 3.817306117981935*^9}, {3.817307189921912*^9, 
   3.8173071903704042`*^9}, {3.817307781609067*^9, 3.8173077822938766`*^9}, {
   3.8173997136294403`*^9, 3.817399723649102*^9}, {3.817399975984412*^9, 
   3.8174000080518312`*^9}, {3.817414176988175*^9, 3.817414177578567*^9}, {
   3.817605365985326*^9, 3.8176053667388163`*^9}, {3.818209885709543*^9, 
   3.818209886669582*^9}, {3.839727321990724*^9, 3.839727332624901*^9}}],

Cell["Example :", "Text",
 CellChangeTimes->{{3.8397049802084637`*^9, 3.839704999360832*^9}}],

Cell["\<\
Suppression factor at 24 degrees for T/T_F=0.11 and E_F/E_R=0.94 :\
\>", "Text",
 CellChangeTimes->{{3.839704874953087*^9, 3.839704973210531*^9}, {
  3.8397050293100224`*^9, 3.8397050441991587`*^9}, {3.839705179247875*^9, 
  3.8397051955931168`*^9}, {3.839705233994288*^9, 3.839705258569729*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"sls", "[", 
  RowBox[{"0.11", ",", "0.94"}], "]"}]], "Input",
 CellChangeTimes->{{3.816444401274846*^9, 3.816444416057328*^9}, {
   3.816708060236553*^9, 3.816708062673461*^9}, {3.816713117973928*^9, 
   3.816713130644359*^9}, {3.816813206935534*^9, 3.8168132072802277`*^9}, {
   3.816827046992669*^9, 3.816827049172257*^9}, {3.816827239987178*^9, 
   3.816827251668414*^9}, {3.816827310272798*^9, 3.816827325479939*^9}, {
   3.81682739238058*^9, 3.8168273989240627`*^9}, {3.816827468563629*^9, 
   3.816827473218226*^9}, {3.816827544688298*^9, 3.816827561012188*^9}, {
   3.816827685029646*^9, 3.816827692394496*^9}, {3.816827800718237*^9, 
   3.8168278081714277`*^9}, {3.8168429070984163`*^9, 3.816842910028714*^9}, {
   3.816842940775674*^9, 3.816842945162648*^9}, {3.816843003361169*^9, 
   3.8168430095665283`*^9}, {3.816843043256557*^9, 3.816843052042832*^9}, {
   3.816843085456503*^9, 3.816843093028359*^9}, {3.816843131989949*^9, 
   3.8168431408020887`*^9}, {3.816843178193919*^9, 3.8168431867747383`*^9}, {
   3.817160316933671*^9, 3.817160325392536*^9}, {3.817169786062645*^9, 
   3.817169790873042*^9}, {3.817306161603039*^9, 3.817306164324144*^9}, {
   3.817306775867951*^9, 3.817306777287994*^9}, 3.8173078536514473`*^9, {
   3.817399752857193*^9, 3.817399759682106*^9}, {3.8174403733375187`*^9, 
   3.817440374188396*^9}, {3.817443082552957*^9, 3.81744314800994*^9}, {
   3.817444296323369*^9, 3.8174443067136927`*^9}, {3.8180506107541523`*^9, 
   3.818050616725951*^9}, {3.8180519533869247`*^9, 3.818051954173135*^9}}],

Cell[BoxData["0.506434635450392`"], "Output",
 CellChangeTimes->{3.839704756304394*^9, 3.8397051428324223`*^9, 
  3.83970526854843*^9, 3.839727949513033*^9, 3.839733326330089*^9, 
  3.839735707923932*^9, 3.8397357836226053`*^9, 3.839735938538096*^9, 
  3.839803731124592*^9}]
}, Open  ]]
},
WindowSize->{808, 755},
WindowMargins->{{19, Automatic}, {Automatic, 31}},
FrontEndVersion->"11.0 for Mac OS X x86 (32-bit, 64-bit Kernel) (July 28, \
2016)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 233, 3, 30, "Text"],
Cell[794, 25, 3107, 85, 180, "Input"],
Cell[3904, 112, 93, 1, 30, "Text"],
Cell[4000, 115, 306, 5, 30, "Text"],
Cell[CellGroupData[{
Cell[4331, 124, 1559, 22, 32, "Input"],
Cell[5893, 148, 275, 4, 32, "Output"]
}, Open  ]]
}
]
*)

